Anisotropic diffusion of electrons in liquid xenon with application to improving 
the sensitivity of direct dark matter searches 



Peter Sorensen 

Lawrence Livennore National Laboratory, 7000 East Ave., Livennore, CA 94550, USA 



Abstract 

Electron diffusion in a liquid xenon time projection chamber has recently been used to infer the z coordinate of a 
particle interaction, from the width of the electron signal. The goal of this technique is to reduce the background event 
rate by discriminating edge events from bulk events. Analyses of dark matter search data which employ it would 
benefit from increased longitudinal electron diffusion. We show that a significant increase is expected if the applied 
electric field is decreased. This observation is trivial to implement but runs contrary to conventional wisdom and 
practice. We also extract a first measurement of the longitudinal diffusion coeflicient, and confirm the expectation that 
electron difilision in liquid xenon is highly anisotropic under typical operating conditions. 
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The liquid xenon time projection chamber has proven 
to be an excellent detector for dark matter searches 
|[T] |2l . The principle of operation is illustrated in Fig. 
[T] and additional details can be found in e.g. ifTSl . In 
standard analyses |l2l|3]|4l|5][6l, the z coordinate of inter- 
action is obtained from the time delay between the elec- 
tron (S2) and primary scintillation (SI) signals. The de- 
lay results from the transit time of the electrons through 
the liquid. The precision is cr^vj < 0.1 cm, where 
(Xe ~ 0.3 fis is the typical width of an S2 pulse, and 
the electron drift velocity is ~ 0.2 cm/jUS [7]. S2 is 
measured via proportional scintillation. The x - y co- 
ordinates of interaction are reconstructed from the hit 
pattern of S2 on the photomultipliers, with precision 
~ 0.3 cm |8|. In this framework, electron diff'usion has 
been safely ignored. 

Recent interest in the possibility that dark matter may 
be OilO) GeV has motivated the XENON 10 collabo- 
ration to explore a new means of analyzing their dark 
matter search data |9|. By ignoring SI and requiring 
only S2 for dark matter candidate events, the detector 
retains full efficiency for nuclear recoils with energy as 
small as ~ 1 keV. This approach is prompted by the fact 
that the detection efficiency for S 1 begins to fall below 
~ 10 keV [10]. In the absence of SI, the z coordinate 
of interaction is indeterminate from the usual method. 
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and incident particle type discrimination (from the ratio 
S2/S1) is lost. However, the width of the S2 signal in- 
creases with drift distance through the liquid, due to the 
diffusion of a swarm of ionized electrons. This allows 
an approximate z coordinate to be assigned to each in- 
teraction, even in the absence of an S 1 signal. Rejection 
of high background event rate edge regions may there- 
fore be partially recovered. Interestingly, such S2-only 
analysis of dark matter search data finds a close parallel 
in the expected sub-keV signals from coherent neutrino- 
nucleus scattering. Experiments which propose ifTTlfTSl 
to search for this signal may also benefit from z coordi- 
nate reconstruction using electron diffusion. 

The rate of transverse electron diffusion in liquid 
xenon is small, about 0.3 cm per meter of drift for 
an applied electric field Ed ^ 1 kV/cm [13]. In this 
regime elastic electron collisions dominate, and solu- 
tion of the Boltzmann equation predicts that diffusion is 
highly anisotropic |14|. We are not aware of any mea- 
surement of the rate of longitudinal (i.e. parallel to Ed) 
electron diffusion. In this brief article we show that the 
longitudinal diffusion coefficient Di of electrons in liq- 
uid xenon may be measured from existing data, and that 
electron diffusion in liquid xenon is highly anisotropic. 
The precision in reconstructed z coordinate using elec- 
tron diffusion is ~ cm, but improves with increasing 
Di. We show that a substantial increase in Di may 
be obtainable with order of magnitude reduction of Ed- 
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Figure 1: Schematic representation of a liquid xenon time projection 
cliamber, showing interaction of a particle X. Primary scintillation 
photons (SI) are generated at the interaction vertex; electrons are also 
generated, and drifted to the gaseous xenon where they create propor- 
tional scintillation (82). A typical event record corresponding to such 
an interaction is shown along the left of the diagram. 



Reducing the uncertainty in reconstructed z coordinate 
would lead to increased efficiency of background rejec- 
tion in S2-only light dark matter searches. The impor- 
tance of this is underscored by the loss of S2/S1 dis- 
crimination in such analyses. 

Consider electrons drifting through Uquid xenon in 
the +z direction under the influence of an applied elec- 
tric field Ed- S2 pulse width cTe from nuclear recoils is 
shown in Fig. |2] as a function of drift time, for events 
in which an SI signal was observed |9|. The nuclear re- 
coil energy range considered is E„r < 30 keV. At these 
energies, the spatial extent of electronic excitation from 
a recoiling xenon nucleus is < 0.1 jjm fTSl, measured 
as nns radius of the track. We therefore treat the initial 
source of no electrons as point-like at f = 0, and write 
their distribution (in co-moving coordinates) at a later 
time t as [16 , 17 1 

,2 



no 



=exp( 



(1) 



V47rDz.f 

We have neglected the transverse dimensions, and note 
that a measurement of Dj is given in fT3l. The mean 
square displacement of the electron distribution is 
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It is more natural to discuss cr, - cr,lvd- The measured 
S2 signal shape is a convolution of the Gaussian elec- 
tron distribution defined by Eq. [T] and the average pro- 
portional scintillation distribution /o from a single elec- 
tron. Proportional scintillation is created in a ~ 0.3 cm 
gap of gaseous xenon under high electric field fW\, with 
a roughly flat probability per unit path length. We model 
/() as a boxcar function with standard deviation ctq. If cr, 
and ctq are similar to within 20%, a criteria that is satis- 
fied for events in which the electrons have drifted for a 
time f > 30 yus, we can write 
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for the measured width of the S2 signal. Equation|4]is 
exact if /o is Gaussian with standard deviation ctq. The 
diff'usion coefficient Di can then be obtained from a fit 
to 
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with (To a free parameter. We find Dl - 12 ± 1 cm^/s 
and o-Q = 0.190 + 0.005 fis for the Ed = 730 V/cm 
XENONIO datiQ fitting only the range t > 40 yKS. This 
is shown in Fig. |2] and the ratio DijDj - 0.15 is in 
good agreement with theoretical expectations IIT4ll . 

The electromagnetic background in liquid xenon dark 
matter detectors is presently dominated by radioactivity 
in the photomultipliers. The primary utility of electron 
diff'usion for light dark matter searches is thus to reject 
events at the top and bottom edges of the active region, 
near the photomultipliers. It can be seen from Fig. |2] 
that in XENONIO, a typical width 0.22 //s is measured 
for events which have drifted 20 jUS. This corresponds to 
a depth of about 3.8 cm below the liquid surface. How- 
ever, the width cr ^ 0.02 //s of the distribution of cr^ 
results in a > 1 cm uncertainty in reconstructed z coor- 
dinate, as shown in Fig. [3] 

In principle, improving the efficiency with which 
edge events are rejected in S2-only analyses could fol- 
low two paths: decreasing the width of the cr^ distri- 
bution, or increasing the electron diffusion. The latter 
appears more accessible. The ratio DilDj approaches 
unity at very low electric field lfT9ll . We therefore ex- 
pect that as Ed is decreased, Di will approach the 
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Figure 2: Width of S2 pulses as a function of drift time through liquid 
xenon. Measurements from XENONIO |9| at fi^ = 730 V/cm and 
Vrf = 0.188 cm//.(s (circles, Icr en'or smaller than data points) and 
XENONIOO |20| at Ej = 530 V/cm and vj = 0.174 cm/ps (squares). 
The total drift distance was 15 cm and 30.6 cm, respectively. D/, 
was obtained from a fit of Eq. |5]in each case. The predicted width 
vs. drift time for two larger values of Z)/. are shown, assuming ctq = 
0. 19 /js and vj = 0. 174 cm/jUS. A simulated distribution of events with 
Gaussian width cr = 0.02 //s are shown for D/. = 40 cm^/s. 



~ 80 cm-/s low-field limit of fT3l. In Fig. |2]we 
show preliminary electron diff'usion data obtained by the 
XENONIOO experiment (stars, uncertainty not shown) 
at Ed = 530 V/cm ll20l . As expected, this data shows an 
increase in D^. 

Dark matter searches with liquid xenon detectors 
have tended to operate at the highest possible value of 
Ed, limited by the maximum voltage that could be ap- 
plied to the cathode. The motivation for this is the ob- 
servation 1 3 1 of improved S2/S 1 discrimination power at 
higher Ed, although other work has found no significant 
improvement ifTSl . For light dark matter searches using 
only the S2 signal, it would be beneficial to operate at 
the lowest practical value of Ed- An optimal value of 
Ed will maximize Di without substantially reducing the 
ionization yield. Studies of the ionization yield of liquid 
xenon for nuclear recoils indicate only a 10% decrease 
if Ed is reduced from 2 kV/cm to 100 V/cm ||2TI . The 
diffusion coefficient at this low drift field is presently 
unknown. 

Suppose Dl were as high as 40 cm^/s at 100 V/cm. A 
simulated distribution of events with width cr - 0.02 /is 
in shown in Fig. |2] In this case, events within ~ 3 cm of 
the edges of the XENONIOO detector would be rejected 
at 3cr by requiring cr^ > 0.31 /is and cTj. < 0.61 fis. 
The z-dependent efficiency of this cut is shown in Fig. 
|4] (a), and results in an effective target mass of 30 kg. 
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Figure 3: Icr uncertainty in z coordinate reconstiucted from S2 pulse 
width varies with Dl as well as z.. Standard z reconstruction from the 
delay time between S2 and S 1 signals results in uncertainty ~ 0. 1 cm. 
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Figure 4: Predicted fraction of accepted events in XENONIOO, with z 
coordinate reconstructed from S2 pulse width, for Di = 16 cm^/s (□) 
and Di = 40 cm^/s (A), cr^ cut bounds for each case are given in the 
text. Dashed lines indicate the accepted range of z coordinate in (2), 
using standard z reconstruction. 



This may be compared with the 40 kg target mass ob- 
tained by the XENONIOO Collaboration from standard 
ts2 - hi z coordinate reconstruction (Fig. |4] dashed). 
Achieving the same level of background reduction with 
Di - 16 cm^/s would require cr^ > 0.27 /is and 
cTj. < 0.39 i-is, and the effective target mass would be 
only 21 kg (Fig.|4) □). 

The sensitivity of liquid xenon detectors to light dark 
matter using the S2-only analysis technique will ulti- 
mately depend on both their ability to reject edge events, 
on the ionization yield per keV of nuclear recoil energy 
f9l and on the S2 signal threshold. The latter is presently 
reported to be about 15 electrons |2 1, significantly above 
the single electron threshold obtained by XENONIO 



||9l. But while XENON 10 was decommissioned in 
2007 with all of its data acquired at Ej - 730 V/cm, 
XENON 100 has the potential to improve an important 
aspect of their sensitivity by turning a single knob. With 
a sufficiently low background event rate in the central 
target, liquid xenon detectors may have the potential for 
unambiguous detection of light dark matter 
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